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Abstract 

The integral of the Wigner function of a quantum mechanical system over a 
region or its boundary in the classical phase plane, is called a quasiprobability 
integral. Unlike a true probability integral, its value may lie outside the interval 
[0, 1]. It is characterized by a corresponding selfadjoint operator, to be called a 
region or contour operator as appropriate, which is determined by the charac- 
teristic function of that region or contour. The spectral problem is studied for 
commuting families of region and contour operators associated with concentric 
disks and circles of given radius a. Their respective eigenvalues are determined 
as functions of a, in terms of the Gauss-Laguerre polynomials. These polyno- 
mials provide a basis of vectors in a Hilbert space carrying the positive discrete 
series representation of the algebra su(l, 1) ~ so(2, 1). The explicit relation 
between the spectra of operators associated with disks and circles with propor- 
tional radii, is given in terms of the discrete variable Meixner polynomials. 
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0. Introduction. 

The problem of bounds on integrals of the Wigner quasidistribution func- 
tion [T] over some subregion or contour of the phase plane for 1-dimensional 
quantum systems, has been introduced and studied as a theoretical issue [2115]) 
with applications in the field of quantum tomography (see e.g. For a 

classical (true) probability distribution, such an integral would always have a 
value in the interval [0, 1], but in the quantum case, the upper and lower bounds 
vary from region to region, and may lie outside the interval [0, 1]. The relevance 
of such studies to quantum tomography, is that these bounds define checks on 
the accuracy of experimental determinations of Wigner functions. It has been 
shown |21 E] that, for a given region or contour, the bounds are determined 
by the maximal and minimal eigenvalues of a corresponding Hermitian opera- 
tor, to be called the region operator or contour operator, respectively, acting 
in the space of square integrable (wave) functions of the quantum systems in 
question. Furthermore, the matrix element of such an operator K between such 
(wave) functions if), gives the value of the non-positive-definite quasiprobability 
integral (QPI), formed by the integral of the Wigner function corresponding 
to ip, over the region or contour that determines the operator. Attention is 
therefore focussed on the spectra of various region and contour operators. 

In what follows, starting from regions and their boundaries in the phase 
plane of a classical mechanical system, we use their respective characteristic 
functions to introduce corresponding region and contour operators. Specifically, 
we introduce the family of commuting operators Kd(o), Kc{a), associated with 
a disk and a circle respectively, with radius a, centered at the origin of the phase 
plane. These operators are diagonal in the Fock space of the number states of 
the quantum harmonic oscillator. Their spectra, which determine the upper and 
lower bounds on the corresponding QPI, show a remarkable richness of structure. 
We show that these spectra, as functions of the radius a, are determined by the 
Gauss-Laguerre (GL) functions 6 , acting as basis functions in vector spaces 
that carry the d = 1/2 positive discrete series representation of the Lie algebra 
so(2,l) ps su(l,l), involving two different differential operator realizations. 
Further, for operators K.D.c{fL a ) that refer to disks/circles with different radii 
£a, £ £ R, we find their eigenvalues are determined by the GL basis functions of 
scaled argument £a, while the latter are expressed as linear combination of the 
GL basis of argument a, with the Meixner discrete polynomials as expansion 
coefficients. 

1. Quantum Mechanical setting. 

A one-dimensional quantum system is described by means of the Heisenberg- 
Weyl (HW) algebra hi, with generators {1, 6, b^} satisfying the canonical com- 
mutation relations bb^ — b^b = [b, 6+] = 1, [b, 1] = = [b>, 1]. In addition to the 
creation operator b\ annihilation operator b and unit operator 1, we also use 
the number operator N = b^b, which satisfies the relations [N, b] — —b, [N, tf] — 
w . An irreducible infinite-dimensional representation v of hi is carried by the 
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so-called Fock space F ~ l2(Z+) = span{e n |n € with 

v(b)e n = y/n&n-u v{b^)e n = Vn + le n+ i, ^(l)e n = e„. 
i/(JV)e„ = ne„, i/(n)e„ = (-l) n e„, (1) 

where we have also introduced the parity operator II for later use. We adopt 
the convention to denote an abstract element and its representative by the 
same symbol. Then, in the representation v, is the hermitian conjugate of 
b, i.e. (ip,bcp)p = (b't/),(p)Fi a s a consequence of the orthonormality relations 
(e m ,e„)j? = $mn- Now we introduce the bounded displacement operators 

D(a) = exp(a*fe- afe 1 ) : C — > B(F). (2) 

Here a is a complex number and a* its complex conjugate. The displacement 
operators define a projective representation of the abelian group of addition in 
the complex plane: 

D(0) = 1, D{ a y = D(-a), (3) 
D(a)D(j3)=D(a + f3)e tax l 3 , (4) 

with a, j3 £ C = T and a x j3 — Re a Im — Im a Re (3. Note that V is endowed 
with the geometric structure of a symplectic space and physically is identified 
with the classical phase space of the quantum mechanical problem dealt with. 
The family of Wigner operators defined as {W'(cc) = D(a)UD(ay\a £ T} pro- 
vides a continously parametrized operator basis for the Hilbert space S(F) of 
Hilbert-Schmidt operators on F due to the relation 

(W(a), W{p)) s = Tr[W{a)^W{f3)} = 5{a - /3) . (5) 

Given A,B£ S(F), so that {A, A) s < oo etc., then 

A = J (A, W{a)) s W(a)d 2 a = J A(a)W(a)d 2 a (6) 

and 

(A,B) S = J^A( a yB(a)d 2 a (7) 

where d 2 a = d(Re ce)d(Im a) /tt. Of special interest are the the density operators 
p £ S(F), which comprise the convex set D of operators satisfying = p. 
p > and Tr[p] = 1. When expanded in the Wigner operator basis, a density 
operator furnishes as expansion coefficient, the corresponding Wigner function : 
W p (a) = (p,W(a)) S - Now Q gives 

Tr(pi) = J A(a)W p (a)d 2 a, (8) 

which is interpreted as the expectation value of the observable represented by 
the operator A, or function A(a), when the system is in the state described 



2 



by the density operator p, or by the Wigner function W p . The set P of 'pure 
state' density operators have the additional projection operator property, i.e 
p 2 = p, and the elements of P are the extremal points of T>, i.e. T> = hull(P). 
For any tp G F, we define ip{x) = ^ n (e n , iji) FH n {x)e~ x / 2 , where H n is the 
(normalized) Hermite polynomial L Q- , so that 

{ih.,ih)F = I Mx)*Mx)dx. (9) 

Jr 

Now if p £ P projects onto ip E F, the Wigner function corresponding to p and 
V> takes the form 

W^p) = ±- [ <p(q + x/2)* x/,(q - x/2) e tpx dx (10) 
27T Jr. 

where q = Re a, p = Ima. The Wigner function was introduced as analo- 
gous to a probability density in an exploration of the extent to which quan- 
tum mechanics can be cast into the form of a true statistical theory ^3 d IH] • 
As is well known and easy to show, W p {q,p) corresponding to some density 
operator p is not in general positive everywhere. It is consequently called a 
quasiprobability function, and some of its previously known bounds, such as 
— (I/ 71 ") < W p (q,p) < (1/tt), as well as some recently found ones [21 El, can be 
used to help quantify its quasiprobability-density character. 

2. Quasiprobability integrals, region operators and contour operators. 

Let Xs ( a ) denote the characteristic function of a given region S in the phase 
space r, and let \ s denote the corresponding operator given by JfjJ), with Xs( a ) 
replacing A(a) in the RHS. We refer to as the region operator corresponding 
to S. If S is compact, is bounded, with a discrete spectrum. 

Considering JSJ) in the case A = x s , we see that the LHS is the expectation 
value of Xs, in the state p, while the RHS is the QPI on S, namely the integral 
of W p over the support of Xs ■ 

It follows that the QPI on S is bounded above and below by the greatest 
and least eigenvalues of Xs , respectively, and also that these bounds are attained 
when p projects onto the corresponding cigenfunctions of Xs- 

A region of special interest is the disk of radius a, centred on the origin, 
with region operator the disk operator x s = Kd(o-). It is known [2] that Krj(a) 
commutes with N for every a. This may be viewed as a consequence of the fact 
that N generates transformations of T which leave every such disk invariant. 
Furthermore, the eigenvalue of K[)(a) on the eigenvector e„ of N has been 
found to be given by means of the Laguerre polynomials L n (x) as, 

\%(a) = 2(-l)" f L n {2x 2 )e- x2 xdx 7 n = 0, 1, 2, . . . . (11) 
Jo 

Generalising the idea of a region operator, we define (3] for every suitably 
smooth contour C in T, a contour operator Xc using jnj), replacing A(a) in 
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the RHS by the generalised characteristic function Xc{a) 01 C- Here Xc( a ) ^ s 
defined by the property that, for every smooth function F(a) on T, 

J F{a)xc(a)d 2 a = J F{a)dl, (12) 

where dl is the element of length along C. We can write 

= / W(a)d 2 a, % = / W(a)dl. (13) 
Js Jc 

In the case that C is the circular boundary of the disk of radius a, centred on 
the origin, we can see that the circle operator = Kc(a) is given by 

K c (a) = lim - (K D (a + e) - K D (a)) . (14) 
e— >o e 

Furthermore, it can now be seen that Kc(a) also commutes with N, and indeed 
with each Ku(b), for every value of a and 6, and that on {e n }„ e z + , Kc(a) has 
the eigenvalue 

\C(a) = ^-\°(a) = 2(-irL n (2a 2 )e- a2 a, n= 0,1,2,.... (15) 



3. Results. 

Our main object now is to clarify the relation between the spectra (|ll|l and 
1151) of commuting concentric circle and disk operators for various radii. 

We begin with some group theoretical preliminaries: Let g denote the Lie 
algebra su(l, 1) with generators {So, S±} in some representation, subject to 
relations 

[S Q , S±] = ±S ± , [S+, S-] = -2S , (16) 

together with Sq — So and S± — S T . The central element C = Sq — 1}S + S- + 
S-S + belongs to the enveloping algebra [7(su(l, 1)). Consider in particular the 
positive discrete series representation T>~^ labelled by a positive integer or half 

integer fc, and carried by a linear space H with formal vectors {d^} n ^z + and 
actions of the operators [5| 

S±d^ = M± 4±i> ^o4 fc) = (k + n)d£\ Cd^ = k(k - l)4 fc ) , (17) 

with = fJ-—(n) = [n(2k + n — l)]i,fi + = H-(n + 1). We will consider two 
different realizations of T>f, called the 7r-realization and cr-realization below. In 
each, the generators are realized by appropriate differential operators acting on 
representation spaces isomorphic to 1 2 {Z + ) . 

7r- Realization. For M real nonzero, the central element and the generators 
of the realization ir : g — ► End{H^) of T>£, are denoted i( M ) = 7r(C), L Q M ^ — 
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tt(So), L± — tt(S±), and act on the Hilbert space ~ L-i ( [0, oo), dr) 



m^{u k i(r)\M G R,k G Z + ,l G {-i,0,i,l,...},r G (0,oo)}. The basis 
vectors are ^U] 

- (-D fc ( J^K) ^ (V5rr)^e-^Lfi(Mr») , (18) 



r(fc + z 



2 V 



with i„(x) the associated Laguerre polynomial 6 . In this representation mod- 
ule the generators are realized by differential operators 

(M) 1 , d 2 1(1 + 1). Mr, 

L± " 4M [ dr 2 + r 2 J + 4 r T 2 { " dr + 2 } ' (2 ° } 

Let ci = i(Z + |), then the basis vectors of H n are eigenvectors of the central 
element with eigenvalues determining the representation i.e. 

L^\r) = d{d-l)u^\r), (21) 
while the action of the remaining generators on the basis vectors reads 

4 M) 4>) = W=«2S ,(r), 4 M) 4>) - (<* + *)4>)- (22) 



a- Realization. Here the central element and the generators of the realization 
a : g — ► End(H a ) of are denoted J = <r(C), Jo = o-(Sq), J± = a(S±), and 
act on the Hilbert space H a ~ L 2 ( [0, oo), r w dr) — spanfe^ (r)\k G R + ,m G 
Z + ,r £ (0, oo)}. The basis vectors are 



e ™ W = 2^^ T ^^exp(-r-)(2r-) fe ^4f- 1 )(2r-), (23) 

with L b a (x) the associated Laguerre polynomial and W — for w > 1. In 
this representation module the generators are realized by differential operators 

Jo = \(w- 2 r 2 - w p 2 r + Zr- W + r w ), (24) 

Ji = \(w- 2 r 2 - w pl - £r- w - r w ), (25) 

J 2 =w- 1 {rp r - l -{w~l)), (26) 

where p r = —i(d/dr + 1/r), £ = ~ 1) — W(W — 1), and the hermitian 
generators J\ — |(J+ — J-) and J 2 — ^(J+ — J-) have been introduced. 

The basis vectors of H a are eigenvectors of the central element with eigen- 
value's determining the representation, i.e. 

Je«(r) = fc(fc-l)e«(r), (27) 
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while the action of the rest of generators on the basis vectors reads 

J± e { n k) (r) = M± 4±iW, Jo 4 fc) (r) = (n + k) (r). (28) 

To make connection with the spectrum of the circle operator Kcifl) we first 
observe that for basis elements of ff w and H a , the relation < f(a),g(a) >h„ = 
< af(a),ag(a) >h„ holds among innner products. The choice of values w = 2 

i.e W — 3/4, leads to the relation «e m 2 4 (a) = wL^Ja), among the basis 
vectors. Then we identify the eigenvalues of Kc(a) with the basis vectors as 
follows: 

A?(o) - aV^el 3) (a) = vW 2) x (a) . (29) 

n, 2 

Proposition 

Let £ G -R+, i/ien £/ie spectra S^' D = {A^ ,z) (^a)|a G R + ,n G <OT<2 S*^' 15 = 

{A„ ,£) (a)|a G R+.n G of the respective circle/disc operators Kc,d{£,o) and 

Kc,d(cl), for concentric circles/disks of radii £a and a, are related by means 
of discrete Meixner polynomials {M n (m,2k;c 2 )\m,n G Z + ,k = 0,^,1,...} as 
follows: 

oo 

A£°(£a) = M» £ c"M„(m, 1; c 2 )A^(a) (30) 
where the coefficients N m , c depend on £. 

Proof. We start by first stating the following, independent of any realization 
of the eigenvalue problem : let X c = — (c + 1/c)Sq + S + + SL , < c < 1, and 
consider £>j: + \ the discrete series representation of g carried by a space with 
basis elements {dn} n ez + - Then 

X c vW = (c- l/c)(k + m)v<£> , meN 7 keR+, (31) 

with 

= £ y^c™M„(m,2fc;c 2 )4 fc ), (32) 

n— 

where {a) n = T(a + n)/r(a) is the Pochhammer symbol and the discrete 
variable Meixner polynomials are defined in terms of the hypergeometric series 
C2|: 

M n {m,2k-c 2 ) = 2 F X [ ~ n 2 ~ m ;l-l/c]. (33) 

[The history of the spectral theory of operators like X c is rich and interesting 
in itself, for it illuminates the relation of special functions with the represen- 
tation theory of su(l, 1), and so provides a mathematical framework for our 
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group theoretical study of region operators. Our ^-realization is based on work 
in ^U]. In J3| > the spectral problem of the operator H = 2J — J+ — J- in 
the positive discrete series representation of su(l, 1) was shown to be related to 
generalized Laguerre polynomials. Following this, the analogous problem was 
addressed f° r H = crJo — J+ — J- , (c S R) , and earlier, for a particular rep- 
resentation, in |15|. The explicit relation to Meixner polynomials in some other 
realizations (cf. our er-realization) was further studied in jllj for the general 
operator H = -(a + 1/c) J - J+ - J_, (0 < a < 1).] 

Since by comparing lfT7|) and (|3*T)l we see that operators So and X c are 
isospectral apart from the factor (c — 1/c), so we expect their formal eigenvectors 
to be related by a unitary transformation. To determine this tranformation we 
first use the Baker-Champell-Hausdorff formula ^B] to rewrite X c in the form 

X c = (c-l/c)e lrS2 S e~ lrS2 , (34) 
where e r = . Its action on the eigenvector vj? then yields 

Soe-^v^J = (k + m)e- trS *v£\ (35) 

which leads to identifying eigenvectors of So and X c up to the normalization 
coefficient N m 

oc 

e**d« =tf m v£> =M m J2<di k \e^d^>di k \ (36) 

71=0 

The right hand side of ijSfil together with yields 

e irS *dW = cfrm Mm £ J^fc n M n (m, 2k; c 2 )d^ , (37) 

V Ti. 

n—0 

where \<j) m \ = 1- 

We first observe that in order to scale the argument of eigenvalues we need 
the dilation operator exp(a^) and that our choice of tt and c-realizations are 
such that their respective operators L 2 = \{L+ — LJ) = — \{ a j^ + \) an d 
J2 = ^{J+ — J-) = \{ a ~ftl + I) are l mear combination of a4-, hence exp(— rL?) 
and exp(irJ2), acting on a function of a, take the respective forms 

cxp(-rL 2 )f(a) = e$/(e*a) and exp(irJ 2 )f(a) = e^/(e'a). (38) 
From this equation and the identification l|29(l . we obtain 

cxp(-r J L 2 + r -)-u^ ] ,(a) = exp(irJ 2 )e^ ) (a), (39) 
2 a ™> 2 

and multiplication of both sides by a^/a yields 

^ (fr) = frv^e^tfa) = (40) 
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Finally, by means of (|36[) and l|37|) and their respective expressions in it, a 
realizations, and with £ = e3 , c = and A/" m = 4>m(^ — c 2 )^c m , where 

specifically for the tt— realization, <f> m = {— l) m , we obtain 

oo 

=M„^c"M„(m,l; C 2 )A^(a). (41) 

From this relation and the spectral decomposition of the circle operator in 
the number state basis, we obtain that 



K c {e^a)= ^ \g(ei a)e m el = ^ (Af m c n M n (m, 1; c 2 )A^(a))e m ej„ , 

in— m— n— 

(42) 

or that 

K c {e^a) = exp(-rL 2 + -)ifc(o) = cxp(ir J 2 )K c (a). (43) 



In view of the simple relation between circle and disk operator spectra, we 
are led to integrate the scaled eigenvalues of (|4*T|l in order to get analoguous 
equations for the eigenvalues of the disk operator, namely 

/>a °° />a 

\g(&)dx=Af m Vc n M n (m,l;c 2 ) / \°(x)dx, 
J ° „ =0 Jo 

OO 

= M„^c"M„(m,l;c 2 )A^(a) (44) 

Below we give a diagrammatic form of the content of the above Proposition: 
the two ladders of spectra for the disc and circle operators with proportional 
radii A„' D (a) and A^ ,£) (£a) respectively, are shown to be related by means 
of coefficents determined by discrete Meixner polynomials, that are given as 
matrix elements of the dilation operator in the Gauss-Laguerre basis functions 
that carry the positive discrete series representation of the algebra su(l, 1) w 
so(2, 1). Also in each spectral ladder the step operators of this algebra J± or L± 
operate to create and annihilate eigenvalues of the respective region or contour 
operators. 
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J_,L_ []J+,L + 



J+,L + T|J_ 5 i_ 



^2' D (£ a ) Meixner pol. 

J-,L- IU+,L + 

\^' D (£a) Meixner pol. 

J-,L- IT J+,L + 



' (£ a ) Meixner pol. 

4. Conclusions. 

The spectra of operators associated with integrals of Wigner functions over 
regions and contours of classical phase space, can have interesting structural 
properties. These operators form a new type of quantum mechanical observ- 
ables, with important mathematical and physical meanings. In the case of discs 
and their boundaries, their eigenvalues, which carry the meaning of the QPI 
formed by Wigner functions over the corresponding discs and circles, behave as 
basis functions of a representation space of a Lie algebra of generators. These 
generators can be used to step up and down the ladder of these eigenvalues. We 
emphasize this unusual feature: the representation is defined on the spectrum of 
eigenvalues, not on the underlying function space on which the region or contour 
operators act. 

Because of their intrinsic interest, as well as their potential importance in 
quantum tomography, it is important to study further the algebraic and func- 
tional properties of such operators. A number of open research problems can be 
mentioned: development of tranformation theory of region and contour opera- 
tors under maps that induce simple geometric transformations on their associ- 
ated region or contour of support; temporal evolution equations and behaviour 
of such operators under Hamiltonian dynamics of their associated quantum 
systems; clarification of their mathematical relation to covariant phase opera- 
tors coming from the problem of quantum mechanical angular operators (see e.g 
|17jh addition rules for region and contour operators associated with many-body 
quantum systems and their relation to geometric manifestations of quantum en- 
tanglement; and finally, the extension of the concept of QPI, region and contour 
operators, and their spectral problems to phase spaces other than the classical 
phase plane, such as the sphere and the continuous and discretized torus. We 
hope to take up some of these problems in future work. 



\C,D, s 



J+,L+ UJ-,L- 
J+,L+UJ-,L. 
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